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Abstract 

We establish the importance of modified Einstein’s gravity (MG) in white dwarfs 
(WDs) for the first time in the literature. We show that MG leads to significantly 
sub- and super-Chandrasekhar limiting mass WDs, depending on a single model 
parameter. However, conventional WDs on approaching Chandrasekhar’s limit are 
expected to trigger type la supernovae (SNela), a key to unravel the evolution¬ 
ary history of the universe. Nevertheless, observations of several peculiar, under¬ 
and over-luminous SNela argue for the limiting mass widely different from Chan¬ 
drasekhar’s limit. Explosions of MG induced sub- and super-Chandrasekhar limiting 
mass WDs explain under- and over-luminous SNela respectively, thus unifying these 
two apparently disjoint sub-classes. Our discovery questions both the global validity 
of Einstein’s gravity and the uniqueness of Chandrasekhar’s limit. 
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Introduction 


The validity of Einstein’s theory of general relativity (GR) has been tested extensively in 
the weak field regime, e.g., through laboratory experiments and solar-system tests. The 
question is, whether GR is the ultimate theory of gravitation, or it requires modification 
in the strong gravity regime. Indeed, it was shown that modified gravity (MG) theories 
reveal significant deviations to the GR solutions for neutron stars (NSs) pQ. As NSs are 
much more compact than white dwarfs (WDs), so far, MG theories have been applied 
only to them in order to test the validity of such theories in the strong field regime. The 
current venture is to show that a MG theory is indispensable in the context of WDs, which 
is a first in the literature to the best of our knowledge. The motivation is the following. 

Type la supernovae (SNela) are believed to result from the violent thermonuclear 
explosion of a carbon-oxygen WD, when its mass approaches the famous Chandrasekhar 
limit of 1.44M 0 , when M 0 is the solar mass. SNIa is used as a standard candle in 
understanding the expansion history of the universe [ 2 ]. 

However, some of these SNela are highly over-luminous, e.g. SN 2003fg, SN 2006gz, 
SN 2007if, SN 2009dc [3j 0], and some others are highly under-luminous, e.g. SN 199lbg, 
SN 1997cn, SN 1998de, SN 1999by [SJ [ 6 j. The luminosity of the former group (super- 
SNela) implies highly super-Chandrasekhar WDs, having mass 2.1 —2.8M 0 , as their most 
plausible progenitors M- While, the latter group (sub-SNela) predicts the progenitor 
mass could be as low as Mq [5j. Moreover, the characteristic lightcurves of these SNela 
are quite peculiar compared to their conventional counterparts. The models attempted 
to explain them so far entail caveats. 

A major concern, however, is that a large array of models is required to explain 
apparently the same phenomena, i.e., triggering of thermonuclear explosions in WDs. It is 
unlikely that nature would seek mutually antagonistic scenarios to exhibit sub- and super- 
SNela. This is where the current work steps in, which unifies the phenomenologically 
disjoint sub-classes of SNela by a single underlying theory. This is achieved by invoking 
a MG theory in WDs. 


Basic equations and modified gravity model 


Let us start with the 4-dimensional action as [7J 


S 


1 

167T 


f(R) + C M 


d A x , 


(i) 


where g is the determinant of the metric tensor g^, Cm the Lagrangian density of the 
matter field, R the scalar curvature defined as R — g' w R iU/1 when R^ v is the Ricci tensor 
and / is an arbitrary function of R (in GR f(R) = R). For the present purpose, we choose 
the Starobinsky model [ 8 j defined as f(R) = R + aR 2 , when a is a constant. However, 
similar effects could also be obtained in other MG theories, e.g. Born-Infeld gravity (e.g. 
0 )- Now, on extremizing the above action one obtains the modified field equation as 


G ^11 + aX uv = 87 rT ; 


fiv 


- fiv'i 


( 2 ) 


where G^ u is Einstein’s field tensor, is the energy-momentum tensor of the matter 
field and X fJU is a function of g^, R^ v and R. 
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Solution procedure 

In order to solve equation (J2J), we adopt the perturbative method (e.g. pH]), such that 
aR « 1. Further, we consider the hydrostatic equilibrium condition: g vr W= 0, 
with zero velocity and V ;i the covariant derivative. Hence, we obtain the differential 
equations for mass M a (r), pressure P Q (r) (or density p Q (r)) and gravitational potential 
0 Q (r), of spherically symmetric WDs, where r is the radial coordinate: the modified 
Tolman-Oppenhcimer-Volkoff (TOV) equations. For a = 0, these equations reduce to 
TOV equations in GR. 

In order to solve the modified TOV equations, we must supplement them with appro¬ 
priate boundary conditions and an equation of state (EoS) obtained by Chandrasekhar 
HU, given by P a = Kp l a +< '^ n \ for extremely low and high densities, where n is the poly¬ 
tropic index and K a dimensional constant. The boundary conditions are M a ( 0) = 0 
and Pa(0) = p c , where p c is the central density of the WD. Note that a particular p c 
corresponds to a particular mass M* and radius R* of WDs. Hence, by varying p c from 
2 x 10 5 grn/cc to 10 11 gm/cc, we construct the mass-radius relation. 


Results 

Figures 1(a) and (b) show that, for a = 0 (GR case), with the increase of p c , M* increases 
and R* decreases, until it reaches a maximum mass M max = 1.4O5M 0 (smaller than the 
Newtonian Chandrasekhar limit of 1.44M 0 ) at p c = 3.5 x 10 10 gm/cc. A further increase 
in p c results in a slight decrease in M*, indicating the onset of an unstable branch. 

Coming to the a < 0 cases, Figure 1(b) shows that for p c > 10 8 gm/cc, the M* — 
p c curves deviate from the GR curve due to MG effects. This reveals that MG has a 
tremendous impact on WDs which so far was completely overlooked. Note that M max 
for all the three cases corresponds to p c = 10 11 gm/cc, an upper-limit chosen to avoid 
possible neutronization. Interestingly, all values of M max are highly super-Chandrasekhar, 
ranging from 1.8 — 2.7M 0 . Thus while the GR effect is very small, MG effect could lead 
to ~ 100% increase in the limiting mass of WDs. The corresponding values of p c are large 
enough to initiate thermonuclear reactions, e.g. they are larger than p c corresponding to 
M max of a = 0 case, whereas the respective core temperatures are expected to be similar. 
This explains the entire range of the observed super-SNela mentioned above Si- 

Table 1 ensures the perturbative validity of the solutions. Recall that we solve the 
modified TOV equations only up to O(a). Since the product aR is first order in a, we 
replace R in it by the zeroth order Ricci scalar R^ 0) = 87r(p^ —3R^), where p ^ and P ® 
are the zeroth order solutions of density and pressure respectively, obtained in GR (when 
a = 0). For the perturbative validity of the entire solution, |a:R( 0 )| max 1 should hold 

true. Next we consider /9tt and g$ / g rr (ratios of g, iu -s in GR and those in MG up 
to O (a)), which should be close to 1 for the validity of perturbative method [13]. Hence, 
|1 — g\^/gtt | max 1 and |1 — gi°J/g rr | ma x •C 1 should both hold true. Table 1 shows 
that all the three measures quantifying perturbative validity are at least 2 — 3 orders of 
magnitude smaller than 1. 

Coming to the a > 0 cases, Figure 1(b) shows that all the three M* — p c curves overlap 
with the a = 0 curve in the low density region. However, with the increase of a, the region 
of overlap recedes to a lower p c . MG effects set in at p c > 10 8 , 4 x 10 7 and 2 x 10 6 gm/cc, 
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Figure 1: Unification diagram for SNela: (a) Mass-radius relations, (b) Variation of p c 
with M*. The numbers adjacent to the various lines denote a/(10 13 cm 2 ). p c , M* and R* 
are in units of 10 6 grn/cc, M 0 and 1000 km, respectively. 


for a = 2 x 10 13 cm 2 , 8 x 10 13 cm 2 and 10 15 cm 2 respectively. For a given a, with the 
increase of p c , M* first increases, reaches a maximum and then decreases, like the a = 0 
case. With the increase of ct, M max decreases and, interestingly, for a = 10 15 cm 2 , it 
is highly sub-Chandrasekhar (O.81M 0 ). In fact, M max for all the chosen a > 0 is sub- 
Chandrasekhar, ranging 1.31 — O.81M 0 . This is a remarkable finding since it establishes 
that even if p c -s for these sub-Chandrasekhar WDs are lower than the conventional value 
at which SNela are usually triggered, an attempt to increase the mass beyond M max with 
increasing p c , for a given a, will lead to a gravitational instability. This presumably will 
be followed by a runaway thermonuclear reaction, provided the core temperature increases 
sufficiently due to collapse. Occurrence of such thermonuclear runway reactions, triggered 
at densities as low as 10 6 gm/cc, has already been demonstrated [hi] . Thus, once M max 
is approached, a SNIa is expected to trigger just like in the a = 0 case, explaining the 
sub-SNela [3 E], like SN 199 lbg mentioned above. Table 2 confirms that the solutions 
for the a > 0 cases are within the perturbative regime. 
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Table 1: Measure of validity of perturbative solutions for a < 0 corresponding to M max 

in Fig. m 


a/(10 13 cm 2 ) 

|aR (0) |max 

11 9tt / 9tt | max 

ll-o (U) /o 1 

I x yrr / yrr \ max 

-i 

0.00184 

0.0016 

0.0052 

-2 

0.00369 

0.0031 

0.0108 

-3.5 

0.00646 

0.0052 

0.0199 


Table 2: Measure of validity of perturbative solutions for a > 0 corresponding to M max 

in Fig. [H 


a/(10 13 cm 2 ) 

|aR (0) |max 

11 9tt / 9tt | max 

\l-a W la 1 

I yrr / yrr \ max 

2 

7.4 x 10 -5 

6.8 x 10“ 5 

2.0 x 10” 4 

8 

7.4 x 10 -5 

6.8 x 10 -5 

2.0 x 10~ 4 

100 

7.4 x 10 -5 

6.9 x 10 -5 

2.0 x 10“ 4 


Conclusions 

Based on a simple /(-R)-model, we show that modifications to GR are indispensable in 
WDs, especially for determining their limiting mass, ft remarkably explains and unifies a 
wide range of observations for which GR is insufficient. We note here that the perturbative 
method is adequate for the present study, as then we have a handle on a characterizing 
our model which cannot be arbitrarily large, allowing it to be constrained directly by 
astrophysical observations [T5j. Hence, depending on the magnitude and sign of a, we not 
only obtain both highly super-Chandrasekhar (for a < 0) and highly sub-Chandrasekhar 
(for a > 0) limiting mass WDs, but we also establish them as progenitors of the peculiar, 
over-luminous and under-luminous SNela, respectively. We further note that even though 
a is assumed to be constant within individual WDs, there is indeed an implicit dependence 
of a on p c , as evident from Figure 1(b), indicating the existence of an underlying chameleon 
effect [16]. Thus, a single underlying theory, inspired by the need to modify Einstein’s 
theory of GR, unifies the two apparently disjoint sub-classes of SNela, which have so far 
hugely puzzled astronomers. The significance of the current work lies in the fact that it 
not only questions the uniqueness of the Chandrasekhar mass-limit for WDs, but it also 
argues for the need of a modified theory of GR to explain the observable universe. 
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